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Abstract 

We derive the equation of state of a dilute atomic Bose gas with an inter- 
atomic interaction that has a negative scattering length and argue that two 
continuous phase transitions, occuring in the gas due to quantum degener- 
acy effects, are preempted by a first-order gas- liquid or gas-solid transition 
depending on the details of the interaction potential. We also discuss the 
consequences of this result for future experiments with magnetically trapped 
spin-polarized atomic gasses such as lithium and cesium. 
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I. INTRODUCTION 



Apart from the intrinsic interest in cooling and trapping atomic gas samples by means 
of electromagnetic fields and their potential use in various high-precision experiments, 
the possibility of achieving Bose-Einstein condensation is also an important motivation for 
studying ultracold atomic gasses in a static magnetic trap. The first steps towards this 
goal were made with atomic hydrogen ||^,|^ and recent progress in lowering the temperature 
of the trapped gas using both conventional and light-induced evaporative cooling 
shows that atomic hydrogen is still a very promising candidate for the ultimate achievement 
of Bose-Einstein condensation. 

Due to their excellent optical properties a substantial experimental effort is presently 
also directed to the alkali- metal vapors cesium [0 and lithium [Q. However, to ensure the 
stability of the condensed phase it is usually (but see below) required that the interaction 
between the atoms is effectively repulsive or more precisely that its scattering length a is 
positive 0. There is no doubt that the above requirement is fulfilled for atomic hydrogen, 
but in the case of atomic lithium the state-of-the-art triplet potential leads to a negative 
scattering length for the doubly spin-polarized |-F = 2,Mp = 2) state ||10[. In contrast, the 
situation for atomic cesium is less straightforward. It has recently been shown that with 
the present knowledge of the interaction potentials it is not possible to determine the sign 
of a for the doubly spin-polarized \F = 4,Mf = 4) state, but that the scattering length 
for the also low- field seeking \F = 3, Mp = —3) state of the lower hyperfine manifold has a 
pronounced resonance structure, which offers the exciting possibility of controlling the sign 
of a by an appropriate choice of the magnetic- field strength |1TT| . 

In view of these circumstances it is of considerable interest to investigate the properties 
of a dilute atomic Bose gas with negative scattering length and to predict what one might 
observe in experiments with magnetically trapped atomic lithium and cesium. In particular, 
it is interesting to see if also for these gasses quantum degeneracy leads to Bose-Einstein 
condensation at sufficiently high densities. However, before we can begin with a detailed 
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discussion of these issues we must be somewhat more precise and reahze that besides the 
scattering length a, an interatomic potential is also characterized by the hard-core radius r^, 
the finite range ry and the well-depth ew- In terms of these quantities the main difference 
between spin-polarized atomic hydrogen and spin-polarized atomic lithium and cesium is 
that in the former case the potential well is so shallow that it admits no bound states, 
whereas in the latter case ew ^ Ti^/mry and many rovibrational states are possible. 

As a result the relationship between the pressure p and the inverse density 1/n at a fixed 
and nonzero temperature T <C eiy/fc^ is qualitatively different and as shown schematically 
in Fig. m At these temperatures atomic lithium and cesium may thus undergo a first- 
order transition from a dilute {nry <^ 1) gaseous phase to a high density {nr^ = 0(1)) 
phase, which is either liquid or solid depending on the details of the interaction potential. 
The coexistence line of these two phases is obtained from a Maxwell construction and is 
therefore such that the areas of region I and II are equal [0. Note that the critical pressure 



Pc{T) found in this manner is always larger than zero, because the area of region I is bounded 
and the area of region II can be made arbitrarily large since for low densities {nA^ <^ 1 where 



A = \j2Txh^ /rakBT is the thermal de Broglie wavelength) the equation of state reduces to 
p = nksT. Hence, for sufficiently low densities the gaseous phase is stable and we are allowed 
to incorporate the infiuence of the interaction by means of perturbation theory around the 
ideal Bose gas. Moreover, at the experimental densities of interest, we are justified in using 
the T-matrix or ladder approximation because nry -C 1 and we only need to include all 
two-body processes in our description of the dilute phase. 

For an adequate determination of the coexistence line, however, a more advanced theory 
is needed which is capable of describing a strongly interacting system at high densities. 
Fortunately, in the case of magnetically trapped atomic gasses we are dealing with a system 
that is expected to have a large free-energy barrier for the formation of a critical bubble of the 
dense phase. Hence, if the density is isothermally increased beyond the point of coexistence, 
a metastable supersaturated vapor will be formed. From an experimental point of view the 
relevant question to be answered is therefore: Can Bose-Einstein condensation also occur 



in the (meta)stable region of the phase diagram (where dp/dn > 0) of an atomic Bose gas 
with negative scattering length? The answer to this question can still be found within the 
framework of the ladder approximation and is the main topic of the present paper. 

In order to do so we must consider the degenerate regime where nA^ = 0(1) and allow 
for the possibility of a second-order phase transition due to Bose degeneracy. Consequently, 
we must first of all identify the appropriate order parameter. This is achieved in Sec. ||. 
Although Sec. |T| is certainly an important part of the paper, the discussion presented there 
is of a somewhat technical nature and makes use of functional methods which are very 
convenient for the derivation of the free-energy density of the order parameters of interest. 
However, after the correct order parameter is found we return in Sec. |T| to canonical methods 
in order to arrive at the final results of the paper in a physically more transparent, and 
perhaps more accessible manner. In particular, we derive in Sec. |T| the equation of state 
of the gas and show that it indeed, in the region of validity of the T-matrix approximation, 
leads to a relationship between the pressure and the inverse density which is qualitatively in 
agreement with Fig. |I|. Most importantly, we are then in a position to answer the question 
mentioned above. We end in Sec. |rV| with some conclusions and with a discussion of the 
experimental relevance of our results. 



We consider a homogeneous system of bosonic atoms in a volume V. Using a functional 
approach to the imaginary time formalism, the Euclidean action of this system is 



where /3 = l/fc^T, fi is the chemical potential, V{x — x') is the interatomic interaction 
potential and possible three-body forces are neglected since we are interested in the dilute 
limit. In addition, the grand canonical partition function is then given by the functional 
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II. ORDER PARAMETER 





integral 

Zg = J d[rW] exp [-]-S[r, ^]) (2) 

over the periodic c- number fields iIj*{x,t) and iIj{x,t). For the ultracold gasses consid- 
ered here the range of the interaction is always much smaller than the thermal de Broglie 
wavelength (ry/A ^ 1), which implies that we can replace the potential V{x — x') by the 
contact interaction Vo6{x — x') with Vq = / dxV{x). Together with the thermodynamic 
potential Q = —{lnZG)/P this completely determines the thermodynamic properties of the 
weakly-interacting Bose gas. 



A. Bose-Einstein condensation 

Nevertheless, to extract useful information from Eqs. (|I|) and (|^) we must resort to 
perturbation theory and thus incorporate the correct symmetry of the homogeneous and at 
high densities metastable state of the gas. Clearly, both objectives can be achieved if we are 
able to identify the correct order parameter. An immediate choice is the thermal average 
{iP{x,t)). However, this choice turns out to be inappropriate for the atomic Bose gas with 
negative scattering length as can already be seen from the following ad absurdum argument: 
Assuming {iIj{x,t)) to be non zero, the application of the usual Bogoliubov theory shows 
that the long- wavelength fluctuations in the order parameter are unstable for a < 0, 
which invalidates the initial assumption. 

A more rigorous argument is based on the free-energy functional for the time and space- 
independent value of the above order parameter 

F[{r),m = vj:7^M)\'\ (3) 

n=l V'-V 

where Fo^"^ is the 2n-point vertex function with all 2n — 1 momentum and frequency ar- 



guments equal to zero |T^,|T3[. Using the T-matrix approximation the latter can easily be 
evaluated in the symmetric (low density or high temperature) phase with the result [|15| 
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T^Q^ = ^S(0; 0)-^ = -fi' , (4a) 
p(4) ^ 2T*^^(0, 0, 0; 2/i') , and (4b) 

rf,^"^ = , for n > 3. (4c) 

In deriving Eq. (|4ED we also made use of the fact that for the momenta hk <^ h/ry and 
Matsubara frequencies Un ^ h/mry of interest the selfenergy hJ]{k; Un) is well approximated 
by the constant 2nT^^(0, 0; 0) = Sirnah^ /m. The selfenergy thus only leads to a shift in the 
chemical potential and we find that fj,' is obtained from the equation of state of the ideal 
Bose gas 

^ = ^E^(^g (5) 

k 

with = fi^l? /2m the kinetic energy of the atoms and A^(e) the Bose distribution l/{e^'^ — l). 
In addition, the many-body T matrix can now be expressed in terms of the two-body T 
matrix T'^^{{], 0; 0) = A'Kah'^/m via ^ 



r^^^(0,0,0;2/iO T2^(0,0;0) e,, - /i' 

Through Eqs. (^) and (^ we observe that for an effectively repulsive interaction 
(a > 0) the free-energy functional has the same form as in the Landau theory of second- 
order phase transitions [0. In particular, the free energy is bounded from below and 
acquires at temperatures below the critical temperature of the ideal Bose gas 



(7) 



a minimum for (ip) = y/noe^^, which determines the condensate density no- However, for an 
effectively attractive interaction (a < 0) this is not true, because the coefficient of the quartic 
term in the free-energy density is negative. We are therefore again led to the conclusion that 
in the case of an atomic Bose gas with negative scattering length {ip{x, r)) is not the correct 
order parameter and that Bose-Einstein condensation in this canonical sence does not take 
place. 



B. Evans-Rashid transition 



Instead, an analogy with the BCS theory of superconductivity suggests that if the 
atoms have an effectively attractive interaction, a phase transition with the order parameter 
A(x, r) = Vo{iP{x,t)iP{x,t)) occurs. To see if this is indeed the case we must derive and 
solve the equation for the equilibrium value of A(x, r) in the ladder approximation, which 
requires the inclusion of fluctuations around the usual BCS (mean-field) theory. As a first 
step towards this goal we perform a Hubbard-Stratonovich transformation by multiplying 
the grand canonical partition Zq by 

l=^fJ d[A*]d[A] exp £^ dr J dx | A(f , r) - Vo^{x, t)^{x, r) ^ j (8) 

and integrating over the fields '?/'*(x, r) and ip{x,T). The latter is conveniently accom- 
plished in Nambu space and requires the introduction of the two-component field 0(x, r) = 
(^/'(x, r), ?/'*(a;, r)) and the corresponding matrix of exact (normal and anomalous) one- 
particle Green's functions G{x, r; x', r') = — (Tr[0(s, T)(f)^{x', r')]), obeying the Gorkov equa- 
tion 



G{x, t; X , t') = h6{x — x)6{t — t') . (9) 

In this manner we arrive at an effective action for the order parameter A(x, r), which is 
formally given by 



S[A\A] = -Tr[lni-G-')]-—J^ dr J dx \A{x,t)\' (10) 

and contains all the desired information on the possibility of a BCS-like phase transition. 
In particular, it can be expanded in powers of A{x, r) by using the Dyson equation = 
Gq^ — E, leading to 

Tr[ln{-G-')] = TrM-G,')] - £ , (n) 

n=l 

and by taking the selfenergy matrix equal to 



^ A(f,r)^ 
yA*(f,r) y 



6{x - x')6{t - t') . (12) 



For our purposes we are especially interested in the quadratic term in this expansion. After 
performing the trace over both coordinate and Nambu space it is found to be 

^(2)[A*,A] = - £^ dr J dx£'^ dr' J dx' A*{x,T)hG^\x,T;x' ,t')A{x' ,t') , (13) 

where the 'noninteracting' Green's function of the order parameter obeys 

Ga^{x, t; X, r) = -^G'o,ii(x, r; x', r')Go,ii(x, r; x', r') + 7^^5(5 - x')6{t - r ) (14) 

or equivalently 



Ga{x, r; x', r') = 2hVoS{x - x')6{t - r') (15) 

dr" / dx" 6*0,11 (a;, t; x", t")Go^u{x, r; x", t")Ga{x" , r"; x', r') . 



Transforming to frequency-momentum space then gives 

1/ g ^ri\Ln + ^f^ ^K/2-k 

having the solution 

Ga{K; ^n) = 2hT^^{6, 0, K; ihVtn + 2/i) . (17) 

In mean-field theory we neglect all fluctuations and as a result conclude that the Evans- 
Rashid transition ||l7i|l9i occurs at a temperature determined by 



/ dr dxGA{x,T;x,T) = GA{0]0) = = 0, (18) 

Jo J , , , ; Av , ; 2hTMB {0,0,0; 2 fi) ^ ' 

which corresponds exactly to the Thouless criterium for the onset of the BCS instability 
p!8| . For a dilute gas, however, mean-field theory is not sufficiently accurate and we must 
also consider the fluctuations. Fortunately, Eq. (|T^ shows how we can apply the T-matrix 
approximation to the Evans-Rashid transition. Indeed, due to this relationship the renor- 
malization of the quadratic term in the action caused by the presence of the | and the | A|^ 
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'interactions', can be accounted for by replacing in the above discussion the noninteracting 
Green's function Go^ii{k] uJn) = h/{ihuJn ~ ej: + fi) by its renormahzed (within the ladder 
approximation) value h/{ihujn ~ ~ ^^{k',uJn) + /i)- Since KL{k\ujn) is well approximated 
by the constant 2nT'^^{0, 0; 0) this implies just a replacement of the chemical potential fi by 
n'. The onset of the BCS-like instability is therefore in the dilute limit determined by 

.1. ^ I. +1e^^^-^""^^-o- (19) 

TA^^(0,0,0;2/) T2S(0,0;0) V ^ - fi' ^ ^ 

Consequently, the critical temperature of the Evans-Rashid transition is slightly above the 



critical temperature of the ideal Bose gas, i.e. Ter = Tq{1 — 0{a/ Aq)) |[T5|| . 

Notice that the above discussion shows that the pair field ip{x, t)iI){x, t) can be used even 
if the interatomic potential has a negative scattering length without having a bound state. 
This is a truly many-body effect which occurs because bosons prefer to scatter into (mo- 
mentum) states that are already occupied. As a result the binding between two particles is 
effectively increased and a resonance, which must always lie just above the continuum tresh- 
old for a to be negative, becomes bound at a density-temperature combination determined 
by Eq. (IID. 



III. STABILITY OF THE GASEOUS PHASES 

In the previous section we have argued that in the degenerate regime the relevant 
order parameter for an atomic Bose gas with negative scattering length is A(x, r) = 
Vo{iP{x,t)iIj{x,t)) . In order to proceed and to discuss whether the corresponding phase 
transition can take place in the (meta) stable region of the phase diagram we must now also 
consider the gas below the critical temperature Ter. This can of course be achieved by the 
functional approach used above, but to make the following more transparent we will from 
now on use canonical methods. 

Denoting the equilibrium value of A{x, r) by Aq and applying the BCS approach in 
combination with the T-matrix approximation, the Hamiltonian of the gas is approximated 
by 
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lA 



^-n2T2^(0,0;0), (20) 



2Vo 

in the Schrodinger picture. After a diagonalisation of this Hamihonian by means of a 
Bogohubov transformation |20| the density n = {ip'^ {x)i/j{x)) and the order parameter Aq = 



Vo{ip{x)tp{x)) are easily calculated, resulting in the equation of state 



(21) 



and the BCS gap equation |T7|JI9 



1 1 1 + 2iv(;^^g) _ 

respectively. 

As will become clear in a moment it is important to note that in the dispersion huj^ 



eg — /x')^ — |AoP of the Bogoliubov quasiparticles there is a minus sign in front of |Aop 
instead of the usual plus sign since we are dealing with paired bosons and not with paired 
fermions. Moreover, the gap equation has an ultraviolet divergence due to the neglect of the 
momentum dependence of the interaction. However, anticipating on the fact that |Ao| is at 
most of O(|nT^-^(0, 0; 0)|) and thus much smaller than fi^/mry, we find from the Lippmann- 
Schwinger equation |^T| for the two-body T matrix that this divergence is cancelled by a 
renormalization of I/Vq to 1/T^^(0, 0; 0). Therefore the gap equation becomes 

1 1 ^ N(huJr) , , 

■ E^-^ = 0, (22) 



T2S(0,0;0) hu^, 

which is free of divergencies. Together with Eq. ( pTD it determines both /x' and |Ao| for a 
fixed density and temperature. 

At high temperatures the gap equation has no solution and we must use | Ao| =0. In that 
case Eq. ( pTD reduces to the expected equation of state for a dilute Bose gas in the normal 
phase (cf. Eq. (|^)). Below the critical temperature Ter, determined by a linearization of 
Eq. (|22|) which correctly leads to Eq. (|19D, the order parameter |Ao| becomes nonzero and 
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the gas is in a superfluid phase with paired atoms. Lowering the temperature, both |Ao| 
and /i' increase, but in such a manner that the gap in the dispersion of the Bogohubov 
quasiparticles decreases since otherwise the system would not be able to accomodate the 
same number of particles. Evidently, this behavior is possible due to the above mentioned 
minus sign in the dispersion relation. 

At a second critical temperature Tbec the gap closes and the number of particles with 
zero momentum diverges, which signals a Bose-Einstein condensation. Below that tempera- 



ture we have |Ao| = — /i' and Tioj^ equals the famous Bogoliubov dispersion \J^\ — '^l^'^k HH 
Furthermore, Eqs. (0) and (|22D are replaced by 

i^^Qy k k ) 

and 

i^-^ + iv^^ = 0, (24) 

determining now yu' and the condensate density tiq for a fixed density and temperature. 
Notice that the mechanism for Bose-Einstein condensation is identical to the mechanism 
causing Bose-Einstein condensation in the ideal Bose gas. In particular, there is no sponta- 
neous breaking of symmetry associated with the second transition and the order parameter 
is no and not {ip{x, r)), which is zero also below Tbec- 

To discuss the question of the mechanical stability of the various phases and consequently 
the experimental observability of the two transitions we must calculate the pressure of the 
gas. Since the Hamiltonian in Eq. ( ]20|) is quadratic this is easily accomplished by evaluating 
the thermodynamic potential = —pV and we find 

J9 = — y |e, - /i' - — - hujr] + + n2T2^(0, 0; 0) 

-^T.^n{l-e-^'-^) , (25) 

k 

where we have again cancelled an ultraviolet divergence in the expression by renormalizing 
1/Vo to 1/T2^(0, 0; 0). The above expression is vahd for all temperatures if we use |Ao| = 
above Ter and |Ao| = — /i' below Tbec- 
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In Fig. we show the behavior of the pressure as a function of inverse density at a fixed 
temperature. At high densities the pressure stays negative, which is unphysical and caused 
by the fact that the T-matrix approximation incorporates the hard core of the interatomic 
potential only in an effective way. For a dilute gas nr^ <^ 1 and this is indeed justified. 
At higher densities, however, the hard core becomes of the utmost importance to avoid 
the complete collapse of the system and must be treated with more care. Nevertheless, we 
can conclude from Fig. ^ that the gas passes in thermal equilibrium through a first-order 
transition from a gaseous phase to a phase with high density. The dense phase cannot be 
described accurately within the framework of the ladder approximation and a more advanced 
theory, capable of describing a strongly interacting system, is required to discuss whether it 
is liquid or solid and also its possible superfluid properties. 

However, as explained in the introduction, for the experimental observability of the 
Evans-Rashid transition we only need to answer the question if it can take place in the 
(meta)stable region of the phase diagram, where dp/dn > 0. It is not difficult to show (by 
comparing Eq. below with Eq. (|I^) that this is never the case. Therefore, we expect 
the following behavior of the atomic system: By increasing the density or lowering the 
temperature, the gas will evolve from a stable to a metastable state until, when dp/dn = 0, 
a point on the spinodal line is reached. Increasing the density further the gas is quenched 
into the unstable region of the phase diagram and will experience a phase separation by 
means of a spinodal decomposition. 



IV. CONCLUSIONS AND DISCUSSION 

In the case of trapped atomic gas samples such a spinodal decomposition should be easily 
visible experimentally because it causes an increase in the density at the center of the trap 
and consequently leads to a sudden increase in the number of two-body relaxation and three- 
body recombination processes that severely limit the lifetime of the system. Therefore, we 
present in Fig. |^ the spinodal line, following from the condition 
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' l^^^^ZiO.o. (26) 



We notice that the degeneracy parameter nA^ is always smaller than C(|) 2.612, implying 
that it is easier to obtain the required condition for phase separation than the condition for 
Bose-Einstein condensation in the case of an atomic gas with positive scattering length. This 
is particularly true for atomic cesium in the \F = 3, Mp = —3) state, which can have large 
negative values of the scattering length by tuning the magnetic field strength. Whether the 
dynamics of the spinodal decomposition is also observable is unclear at this moment because 
it requires a detailed study of the growth of the liquid or solid phase taking into account 
the release of latent heat and the enhanced importance of inelastic collision processes that 
lead to decay of the sample. 

Summarizing, we have studied the dilute atomic Bose gas with negative scattering length 
and have shown that although Bose-Einstein condensation can in principle also occur in this 
case, it does not take place in the gas phase. In contrast, we argue that the gas separates into 
a normal gas and a (possibly superfluid) liquid or sohd. For magnetically trapped atomic 
lithium or cesium the most likely scenario is phase separation into a gas and a solid since 
both the thermal energy as well as the energy due to zero-point motion are much smaller 
than the depth of the potential well. However, a more elaborate discussion is necessary 
to confirm this conjecture. In any case the separation proceeds by means of a spinodal 
decomposition and gives a clear experimental signal if one monitors the decay of the atomic 
density. 
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FIGURES 

FIG. 1. Qualitative picture of the relationship between pressure and inverse density at a fixed 

and sufficiently low temperature for 1) atomic lithium or cesium and 2) atomic hydrogen. Note 
that possible cusps in the curves due to symmetry-breaking phase transitions are not shown here. 

FIG. 2. Pressure as a function of inverse density for a fixed temperature such that |a/A| = 10~^. 
The inset shows the cusp in the pressure caused by the Evans-Rashid transition. 

FIG. 3. Spinodal line in the (nA^, na^) plane, separating the (meta)stable and unstable regions 

of the phase diagram. The dashed line corresponds to the condition for Bose-Einstein condensation 
in a Bose gas with repulsive interactions. 
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